Dynamics of Lieb-Liniger gases 
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It is proved that the Lieb-Liniger (LL) cusp condition implementing the delta function interaction 
in one-dimensional Bose gases is dynamically conserved under phase imprinting by pulses of arbitrary 
spatial form and the subsequent many-body dynamics in the thermodynamic limit is expressed 
approximately in terms of solutions of the time-dependent single-particle Schrodinger equation for a 
set of time-dependent orbitals evolving from an initial LL-Fermi sea. As an illustrative application, 
generation of gray solitons in a LL gas on a ring by a phase imprinting pulse is studied. 

PACS numbers: 03.75.-b,05.30.Jp 



Recent advances and experiments on ultracold atomic 
vwors in atom waveguides I", "3, II II IE S II [s*, 9*, To , Tl', 

I12L 1 131 and potential applicability to atom intcrfcrome- 
try |l4| and integrated atom optics 0, 0| create a need 
for accurate theoretical modelling of such systems, and 
this is underscored by recent experimental achievement of 
optical atom interferometric structures of Mach-Zchnder 
and Michelson type At sufficiently low tempera- 

tures and densities and high transverse frequencies loq 
the transverse degrees of freedom of a Bose-Einstein con- 
densate (BEC) in an atom waveguide are frozen in the 
ground transverse mode and the dynamics is described 
by an effective ID Hamiltonian with delta function inter- 
actions 0, 01 . In the spatially uniform case the exact 
gro und state was found in 1963 by Lieb and Liniger (LL) 
[ijj- The properties of this LL gas are determined by a 
dimensionless parameter 7 — mguj / nh^ where n is the 
ID density and the effective ID coupling constant 
is related to the 3D s-wave scattering length a by 
giD = 2ah^/mtl[l - (C/\/2)(a/4)] with 4 = ^/h/muj^ 
the transverse oscillator length and C = 1.4603 .... If 
7^1 the dynamics reduce to those of a I D g as of 
hard core, or impenetrabl e, p oint bosons 0, Il8l |. the 
"TG gas" ,21, 23,, .23, iM Hi- This is a model for 
which exact many-body energy eigensolutions and dy- 
namics were found using a mapping from the Hilbert 
space of eigenstates of an ideal gas of spinless fermions 
to that of mariy-body ei gens tates of hard core bosons 
lilillMliSISillllllli. The to regime is only a 
tiny part of the LL regime, most of the latter having al- 
ready been reached experimentally 0, 0, 0, 0, 0| , and 
it is desirable to extend the techniques applicable in the 
TG regime to the whole LL regime so as to reduce the 
many-body dynamics of free and trapped LL gases to so- 
lution of the time-dependent single-particle Schrodinger 
equation for a set of time-dependent orbitals evolving 
from initial orbitals comprising a LL-Fermi sea. It will 
be shown herein that although this cannot be done ex- 
actly for all values of 7 and all time intervals, it can be 
done in good approximation for limited time intervals for 
all 7, and exactly in some special cases. 

LL cusp condition and energy eigenstates: The A^- 



particle energy eigenstates satisfy HoijjBa = EaijjBa 
where V'Bq is symmetric (Bose). LL considered a uni- 
form system and used dimensionless units such that 
h = 2m = 1, yielding a Hamiltonian 
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where 2c = gio 
conditions ipBai' 



2^n. They used periodic boundary 



■ Xj ± L ■ ■ ■) 



■). The A^- 



particle configuration space decomposes into A^! permu- 
tation sectors, each defined by one particular ordering 
of numerical values of the xi, • • • ,xn- In view of the 
Bose symmetry of the ijjBa, it is sufficient to define them 
only in the fundamental permutation sector wherein 
xi < X2 < ■ ■ ■ < xn, then extending to all A! sectors 
by symmetry. The cusp condition following from the 
S{xj — Xk) interactions is |l9j 
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dx, dx, ' ^^"1---'= + 



cV'BaU,-=a:fc+ • (2) 



Define 13 j — d/dxj+i — d/dxj — cl where 1 is the unit 
operator. Then in £Hi 

Bj1pBa\x, + i=x,+ = ■ (3) 

In D\i the energy eigenstates satisfying this are 

• • • , :r^) = ^2= J2 ap{K^)e^^P^-y^ (4) 

and their energies are Ea — ^J- Here is an N- 
dimensional vector (fci, A:2, ■ • • , fc^r) with fci < A:2 < • • ■ < 
kjsi, the sum runs over the A^! permutations P of this 
set, and X is the A^-dimensional position vector X = 
(xi, X2, • • • , xat), assuming xi < X2 < • • • < xn- There 
are infinitely many such choices /ci, • • • , fcjv, and a labels 
each particular choice. The coefficients are |3 1^ 



a (K ] = e TT in + ijPki- Pkj) 



(5) 
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where Pkj is the image of kj under P and ep is the parity 
±1 of P. This differs from the ap of LL Eqs. (2.12) ff. 
only by a constant phase. 

Ahhough the symmetrical extension of i^Ba to all per- 
mutation sectors is bosonic (totally symmetric under 
permutations of xi, ■ ■ ■ ,xjv), it can be shown that 
it is nevertheless aniii symmetric under permutations of 
ki^ - ■ ■ ,kM- Thus V'-Ba must be built from N singly- 
occupied orbitals. The ground state V'so is obtained by a 
particular choice Kq consisting of the smallest N differ- 
ent kj consistent with the periodic boundary condition 
on the "iliBa * "LL- Fermi sea". 

Thermodynamic limit: The allowed kj were deter- 
mined in the LL paper by requiring that ipsaixi, • • • , xn) 
be periodic in each Xj with periodicity length L, lead- 
ing to {kj+i — kj)L = J27=ii^sj — ^sj+i) + ^TTUj where 
6sj = — 2 tan""'^ [(fcs — kj)/jn] and the Uj are any inte- 
gers > 1, whose values distinguish between all the TV- 
particle energy eigenstates 19, 35]. The ground state is 
given by Uj = 1 for all j. It follows that the individ- 
ual orbitals e**"'^^ in Q are not L-periodic; the kj are 
not actual particle momenta but rather nonuniformly 
spaced quasimomenta [l^. Since our goal is reduction 
of the many-particle time-dependent Schrodinger equa- 
tion (MBTDSE) to the single-particle time-dependent 
Schrodinger equation (SPTDSE), this leads to difficulties 
related to nonhermiticity of the single-particle Hamilto- 
nian and an ill-posed initial value problem for the orbital 
evolutes (j)j{x,t) on the finite interval —^<x<^. 
These can be mitigated in the thermodynamic limit 
A'' oo, L oo, N/L n ~ const, wherein both the 
quasimomenta kj and the actual particle momenta k^ 
become dense and every quasimomentum can be approx- 
imated with arbitrary accuracy by the k^ — v2'!i/L with 
V ~ 0, ±1, ±2, • • •, appropriate to periodic boundary con- 
ditions. The initial orbitals e^^"'^ j^fL and their evolutes 
(l)v{x,t) are then orthonormal on < x < ^. In this 
limit any thermodynamically intensive or extensive quan- 
tity expressible as a sum of contributions of "quasior- 
bitals" evolving from the quasimomentum plane waves 
^ikjx approaches a similar expression expressed in terms 
of contributions of L-periodic orbitals evolving from L- 
periodic plane waves e^^"^ . In doing this it is necessary to 
take into account the nontrivial density of states of the 
allowed quasimomenta. For example, the ground state 
energy in the thermodynamic limit is '}2!^J'^~vq P^^'^)^t 
where (i{ku) is the ratio of the LL ground state den- 
sity of states to its limiting value l/27r in the TG limit, 
and z/Q = LkLLp/2iT where khLF, the wave vector at the 
top of the LL-Fermi sea, is determined by the condition 
j':_llU{k)dk = 2nn^. 

Response to phase imprinting: Solitons have been gen- 
erated in experiments on cigar-shaped BECs by appli- 
cation of phase- imprinting pulses |36l Is^ . We have re- 
cently studied 28] the exact dynamics of such a process 



in the TG limit, as well as phase-imprinting pulse gen- 
erated interference effects in a TG gas on a ring jSOl IS ij 
and in a double-X atom interferometer |32l ]. A similar 
LL model will be studied here. Consider a LL gas on 
a ring of circumference L initially in its ground state, 
and represent the positions of the A^ atoms by a ID co- 
ordinate — -J < X < -J measured around the circumfer- 
ence. Suppose that this system is subjected at i = 
to a phase-imprinting pulse v{x,t) — — (5(i)S'(a;). Denot- 
ing the orbitals of the LL solution by e'^^^^, the orbitals 
just after the pulse are [s^ <l)j{x,t = 0+) = e*'S'(^)e*''j^, 
and by (0J just after the pulse iPb{xi, ■ ■ ■ ,XN;t = 
0+) = ipBaixi, - ■ ■ ,XN)exp{iY,^=iS{xj)). Supposing 

that xi < a:2 < • • • < xn, and using the Bj of Eqs. 
and (j2J), one finds 

Bj^BiO+) - i^B{0+)[S'ix,+i)~S'{x,)] + e'^B,i,Bc. 

(6) 

with S = '^jS{xj). This vanishes as Xj+i — > Xj-\-. It 
follows that the cusp condition need not be imposed as a 
constraint during the pulse, being exactly valid at i = 0+ 
independently of the spatial form and strength of the 
pulse. 

Dynamics: Consider next free propagation according 
to Ho following such a pulse. This is generated in by 
the kinetic energy alone, i.e., by Eq. ^ with the delta 
function interactions omitted. Let if^Bixi, ■ ■ ■ ,xj^;t) be 
a wave function differing from only by replacement of 
the orbitals e'^^J^ by the solutions 4>j{x,t) of the single- 
particle TDSE [id/dt-\- d^/dx^](j>j{x,t) = reducing to 
(f>j{x, t = 0-I-) as i ^ 0-1- . Then ipBi^i, ■ ■ ■ ,XN;t) satis- 
fies the MBTDSE [ihd/dt- Ho{t)]ijB = in the interior 
of fHi for f > and reduces to the previously discussed 
wave function ij]B{xi, ■ ■ ■ , xn; t — 0-f ) as t ^ 0+. How- 
ever, a subtle complication enters here: Although the 
derivative operators in Bj of Eq. ^ commute with the 
kinetic energy, the ratio of third to second derivative of 
the wave function is different from the ratio of first to 
zeroth derivative required by the cusp condition [l^l23|. 
so a wave function propagated as above does not in gen- 
eral satisfy the cusp condition for t > 0. However, it is 
satisfied trivially for 7^0, and also for 7 ^ 00 (TG 
limit). Furthermore, by temporal continuity of the solu- 
tion of the TDSE it is satisfied in good approximation 
for a nonzero range of times t > for all 7, assuming as 
in © that V'b(- ■ • ; * = 0) is an energy eigenstate 
By a temporal phase mixing argument the length of this 
range should increase with decreasing values of AE/E, 
suggesting that constraint-free propagation as above is a 
reasonable approximation for AE/E ^ 1 where E is the 
energy and AE is its dispersion. 

Single particle density and soliton propagation: Sup- 
pose now that the initial A^-particle wave function 
iPb{- ■ ■ ',t = 0+) differs from the LL ground state i/'so 
at most by multiplication of the LL orbitals e^^i^ j^fh 
by e*'^^^-' as in the previously discussed case of phase im- 
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printing. Then ipsi' ■ ■ \ t = 0+) will satisfy the LL cusp 
condition exactly. Suppose furthermore that /^E/E ^ 1 
where A£' is the energy uncertainty of • • ; ^ = 0+). 
Then to a good approximation one has in IHi 

tpBixi, ■ ■ ■,XN;t) 

= -7=^ap{Ko)(t>Pi{xi,t) ■ ■ ■ 4)pn{xn :t) (7) 
V jV. p 



where the (j)j{x,t) are evolved by the single-particle 
TDSE starting from 0j(a;,O+) = e^'^^'^e''^'-^^ /y/Z. 
The single particle density at time t is n{x,t) — 
N J • ■ ■ J \ipBix, X2, - ■ ■ , xn] t)\'^dx2 ■ ■ ■ dxN where the in- 
tegrals run over one periodicity cell —L/2 < xj < L/2 
and the thermodynamic limit is assumed so that pe- 
riodic exponential orbitals are used. Upon substitu- 
tion from Q one obtains a double permutation sum 
J2pp'- Consider first the terms P ^ P' . Recall that 
holds only in 5Hi and is to be extended to other sec- 
tors by Bose symmetry. Then since P' differs from P 
by at least one interchange, and in most cases by many, 
integrals of the form j_^j^^^(j)*{xm,i)(j)i{xm,t)dxm and 

J^^^'^ (j>*{xm,t)(j)e{xm,t)dxm arise, where j, I, m, and p 
are all different. The phases of these integrals have de- 
nominators L following from the normalization of the 
4>j and their phases oscillate as a function j — £, 
and this phase variation is compounded by oscillation 
of the phases of the factors a*p{KQ)apr{KQ) from Eq. 

One therefore expects that in the thermodynamic 
limit the net off-diagonal contribution will be small com- 
pared with the diagonal ones and will be dropped, an 
approximation in the same spirit as the random phase 
approximation. The diagonal contributions P = P' be- 
have very differently. Since (0 holds only in 9li and 
is to be extended to other permutation sectors by Bose 
symmetry, the integrations over X2, ■ ■ ■ , xn have different 
limits for each different ordering of x, a;2, ■ ■ • , xn- How- 
ever, summing over all permutations, collecting coeffi- 
cients of and performing appropriate permuta- 
tions of the names of integration variables, one finds that 
the integrations over X2, - ■ ■ ,xn sum to a multiple of 
" I-l/2^^n\(I>2{x2)\'^ ■ ■■\(I}n{xn)\'^ = 1, and 
the same holds for the coefficients of the other 
Noting also that |ap(ifo)P = 1 by jSJ and that there 
are {N — 1)1 permutations P such that PI = j for each 
j = 1, ■ ■ ■ , N, one finds that the final expression collapses 
simply to the sum of absolute squares of all N orbitals. 
Recalling that approximation of the nonperiodic LL ex- 
ponentials e^'^J^ by the periodic exponentials e**^"^ in the 
thermodynamic limit requires inclusion of a density of 
states factor, one obtains 

n{x,t)^ J2 0iK)\Mx,t)\^ (8) 

where and f3{k^) were defined in the previous discus- 
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FIG. 1: Scaled density n{x,t)L versus scaled position x/L 
around the ring using 65 orbitals, at time t after a n phase 
imprint imposed for —L/2 < x < 0. The middle half — 1//4 < 
X < L/4 of the periodicity cell is shown. Line 1: 7 = 1000 
(TG limit), L = 1, t = .OOlr. Line 2: 7 = 1, L = 2.19, 
t = .005t. Line 3: 7 = .05 (Bogoliubov regime), L = 7.73, 
t = .06r. Here r = l/27r in the units used (fi — 2m — 1, mean 
density n — 65). 

sion of the thermodynamic limit. 

To apply this to gray soliton generation consider a ring 
of circumference L with circumferential coordinate x with 
—L/2 < X < L/2 and impose a phase shift tt on the left 
half —L/2 < X < 0. The density of states was determined 
by solution of Eqs. (3.18) and (3.20) of LL iiaiSS], and 
the TDSE propagation was effected by EFT of the ini- 
tial orbitals ^^f^^^^^six) ^f^gj. -^Jiich the resultant Fourier 
components have trivial time dependence e~ . The 
density functions showing soliton propagation are shown 
in Fig. ^for the cases 7 = 00 (TG limit), 7 = 1 (inter- 
mediate regime), and 7 = 0.05 (Bogoliubov regime) for 
vq — 32, corresponding to a LL-Fermi sea of 65 orbitals. 
The periodicity cell lengths L were adjusted to keep the 
mean density (65 in the units used) and hence the LL- 
Fermi momentum fixed, using Eqs. (3.14)-(3.21) of LL 
[1^ . The shape of the density profile is quite insensitive 
to 7 apart from the indicated scaling of the periodic- 
ity length L and observation time i, in spite of strong 
peaking of the density of states about low momenta as 
7 varies from the TG regime 7 ^ 1 to the Bogoliubov 
regime 7^1. This insensitivity is consistent with scal- 
ing of the allowed momenta ku = 2-k / L inversely with L 
and scaling of observation times t in their propagation 
factors e^*'^"'* directly with L^. 

One can define a LL-Fermi sea quasimomentum distri- 
bution function f{k) as the absolute square of the Fourier 
transform of each orbital, weighted by the LL quasimo- 
mentum density of states and summed over all orbitals. 
This was evaluated with the aid of an EFT code and is 
shown in Fig. |21for the case 7 = 1 before and just after 
the TT phase imprint. The perturbation promotes particles 
just under the LL-Fermi surface to states just above the 
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FIG. 2: Quasimomentum distributions f{k) for j — 1, nor- 
malized to total number of particles nL, with n and L the 
same as in Fig. 1. Line 1: Unperturbed LL-Fermi sea. Line 
2: Perturbed sea after a vr phase imprint applied to half of 
the ring. 

surface, indicating that Pauli blocking is still quite effec- 
tive at 7 = 1, to be expected since the density of states 
is still almost flat. The quasimomentum distribution is 
not the same as the true momentum distribution, which 
is peaked about fc = even in the TG hmit UJ and 
becomes increasingly peaked as 7 decreases |2^. Never- 
theless, the energy is given exactly by = 
both for the LL ground state 19] and the LL excited 
states (3S] . It is thus reasonable to approximate the en- 
ergy of the phase-imprinted state by the same simple ex- 
pression. Furthermore, since one starts from the ground 
state which has zero dispersion, the energy dispersion af- 
ter the pulse should be of the same order as the mean 
excitation energy Eafter — Ebefore- The FFT evaluation 
yields {Eafter - Ebef ore) /Ebe fore = 0.083 < 1, Suggest- 
ing that the soliton density profiles of Fig. ^ are good 
approximations to exact profiles which are not currently 
calculable. 

Discussion: The LL solitons exhibit density ripples ab- 
sent from Gross-Pitaevskii solitons even for 7 as small as 
.05, being many-orbital features stemming from fermion- 
ization which is present for any nonzero 7, no matter 
how small [l9l |3J| . Even the ground state of a trapped 
LL gas is currently unknown, so it is not clear whether or 
not such ripples will be present in soliton-likc dynamics 
of a trapped LL-gas, an important topic for future in- 
vestigations. The approach used for 6{t) pulses could be 
extended to arbitrary potentials v{x,t) by split-operator 
propagation of orbitals. Could it be extended to the 
trapped LL gas? 
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